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Quantum interference arising from superposition of states is a striking evidence of the validity of Quantum 
Mechanics, confirmed in many experiments and also exploited in applications. However, as for any scientific 
theory. Quantum Mechanics is valid within the limits in which it has been experimentally verified. In order to 
extend such limits, it is necessary to observe quantum interference in unexplored conditions such as moving 
terminals at large distance in Space. Here we experimentally demonstrate single photon interference at a ground 
station due to the coherent superposition of two temporal modes reflected by a rapidly moving satellite thousand 
kilometers away. The relative speed of the satellite induces a varying modulation in the interference pattern. 
The measurement of the satellite distance in real time by laser ranging allowed us to precisely predict the 
instantaneous value of the interference phase. We then observed the interference patterns with visibility up to 
67% with three different satellites and with path length up to 5000 km. Our results attest the viability of photon 
temporal modes for fundamental tests of Physics and Quantum Communications in Space. 


Introduction - Quantum interference has played a crucial 
role to highlight the essence of Quantum Mechanics since the 
Einstein-Bohr dialogues at the end of the Twenties Q. In¬ 
deed, it originates when alternative possibilities in a quan¬ 
tum process are indistinguishable, like in the case of indi¬ 
vidual particles that may be simultaneously in more than one 
place as in the well-known Young double-slit experiment fT]. 
Quantum interference has been observed with photons ||2l[3l. 
but also with electrons ||4|, neutrons jS] and even with large 
molecules with masses exceeding 10000 amu 0. One of the 
main challenges in Quantum Physics is establishing if funda¬ 
mental bounds to interference exist: for instance, can quantum 
interference be measured by observers in relative motion and 
at arbitrary large distance? To our knowledge, the longest path 
on which interference at the single photon level was tested is a 
307 km fiber link on ground with fixed sender and receiver Q. 
Classical interference in free-space has been observed in grav¬ 
itational wave detectors based on a Michelson interferometer 
with 4 km long arms and using a laser beam with kilowatt 
power 0. 

Here we demonstrate interference at the single photon level 
along satellite-ground channels by exploiting temporal modes 
of single photons. To this purpose, we exploited a coher¬ 
ent superposition between two single-photon wavepackets on 
ground and observed their interference after the reflection by 
a rapidly moving satellite at very large distance with a total 
path length up to 5000 km. The varying relative velocity of 
the satellite with respect to the ground introduces a modu¬ 
lation in the interference pattern which can be predicted by 
special relativistic calculations, as explained below. 

Description of the experiment - In our scheme, a coher¬ 
ent state I flyout) in two temporal modes is generated at the 
ground station with an unbalanced Mach-Zehnder interferom¬ 
eter (MZI), sketched in the lower panel of Fig. The delay 
At ~ 3.4 ns between the wavepackets of the two modes 
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FIG. 1. Scheme of the experiment and satellite radial velocity. In 

the top panel we show the measured radial velocity of the Beacon-C 
satellite ranging from —6 km/s to -|-6 km/s as a function of time dur¬ 
ing a single passage. In the bottom panel we show the unbalanced 
MZI with the two 4/-systems used for the generation of the state 
and the measurement of the interference. Light and dark green lines 
respectively represent the beams outgoing to and ingoing from the 
telescope. In the inset, we show the expected detection pattern: the 
number of counts Nc in the central peak varies according to the kine¬ 
matic phase ip imposed by the satellite. Right photo shows MLRO 
with the laser ranging beam and the Beacon-C satellite (not to scale). 
The phase ip{t) depends on the satellite radial velocity as described 
in the text. 


corresponds to a length difference between the two arms of 
I = cAf ~ 1 m (c is the speed of light in vacuum) and it 
is much longer than the coherence time ~ 83 ps of each 
wavepacket (we used the convention of 191 for the definition of 
Tf,). Using a telescope, the state |4'out) is directed to a retrore- 
flector placed on a satellite in orbit. The satellite retroreflec- 
tors redirect the beam back to the ground station, where it is 
collected and injected into the same MZI used in the uplink. 
After the reflection by the satellite and the downlink attenu- 
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ation, the state collected by the telescope can be written as 
= (l/'\/2)(|S') — namely as a a superposition 

of two single-photon wavepackets [S') and \L) (the quantum 
state written above corresponds to the re-normalized single¬ 
photon part of the state received at the telescope). We note that 
the above superposition is also known as time-bin encoding, 
and it is used for fundamental tests of Quantum Mechanics 
UMia, for Quantum Information applications such as quan¬ 
tum key distribution (QKD) lfT4lfT5l along optical fibers ifThl - 
[TS] and for increasing the dimension of the Hilbert space in 
which information can be encoded 1191 . 

As we now explain, the relative phase ip{t) is determined 
by the satellite instantaneous radial velocity with respect to 
ground, Vr{t). Indeed, at a given instant t, the satellite motion 
determines a shift 6r{t) of the reflector radial position, during 
the separation At between the two wavepackets. This shift 
can be estimated at the first order as Sr{t) « Vr{t)At, and 
its value may reach a few tens of micrometers for the satel¬ 
lites here used. For instance, in the top panel of Fig. [T] we 
show the value of Vr{t), that ranges from -6 to 6 km/s for 
the selected passage of the Beacon-C satellite. Therefore, the 
satellite motion imposes during reflection the additional kine¬ 
matic phase (p{t) « 2(5r(f)(27r/A) between the wavepackets 
\L) and jS”), where A is the pulse wavelength in vacuum (see 

Fig- Hi- 

A single MZI for state generation and detection intrinsi¬ 
cally ensures the same unbalance of the arms and avoids ac¬ 
tive stabilization, necessary otherwise with two independent 
interferometers. As detailed in Appendix [A] two 4/-systems 
realizing an optical relay equal to the arm length difference 
were placed in the long arm of the MZI. The relay is re¬ 
quired to match the interfering beam wavefronts that are dis¬ 
torted by the passage through atmospheric turbulence: other¬ 
wise, the latter may cause distinguishability between the two 
paths, washing out the interference. The MZI at the receiver 
is able to reveal the interference between the two returning 
wavepackets. At the MZI outputs we expect detection times 
that follow the well known three-peak profile (see Fig. [^l: the 
first peak represents the pulse [S') taking again the short arm, 
while the third represents the delayed pulse \L) taking again 
the long arm. In the central peak we expect indistinguishably 
between two alternative possibilities: the \S) pulse taking the 
long arm and the \L) pulse taking the short arm in the path 
along the MZI toward the detector. The signature of inter¬ 
ference at the single photon level is then obtained when the 
counts in the central peak differ from the sum of the counts 
registered in the lateral peaks. 

To measure the interference we used a single photon detec¬ 
tor (PMT) placed at the available port of the MZI, as shown in 
Fig. 1. For a moving retroreflector, as detailed in Appendix^ 
special relativistic calculations show that the probability of 
detecting the photon in the central peak is given by 


Pcit) = ^[1-V{t) cos (p{t)] , 


( 1 ) 
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FIG. 2. Kinematic phase and interference pattern. Top panel: 
we show the measured satellite distance and the predicted kinematic 
phase estimated by eq. 0 as a function of time for a passage of 
the Ajisai satellite. Shaded area represents the temporal window of 
data acquisition. Bottom panel: kinematic phase ip{t) and theoreti¬ 
cal probability Pc{t) in the shaded area. The interference pattern is 
modulated according to the value of ip{t) determined by the satellite 
velocity. 


with 


ip{t) 


2/)(f) 27rc ^^ 

1+m A 


- ^ V. (*) r, r At /3(f) 

V(f) = e = e i+ 3 (t)) 


( 2 ) 


We note that for a retroreflector at rest we expect Pc = 0. 
The parameter /3(f) is defined as /3(f) = The above 

relation is obtained by time-of-flight calculations together 
with the Doppler effect that changes the angular frequency 
of the reflected pulses from ujq = ^ to We note 

that the first order approximation of eq. ([^ gives the phase 
(/3(f) « 47rUr(f)Af/A above described. The theoretical vis¬ 
ibility V(f) is approximately 1 since the /3 factor is upper 
bounded by 3 • 10“® in all the experimental studied cases, 
while the ratio Af /tc is of the order of 10^. 

Experimental results - We realized our experiment at the 
Matera Laser Ranging Observatory (MLRO) of the Italian 
Space Agency in Matera, Italy, that is equipped with a 1.5 
m telescope designed for precise satellite tracking and which 
acted as ground quantum-hub for the first demonstrations of 
Space Quantum Communication (QC) MM- The pulses 
used to prepare I'I'out) are generated by a mode-locking laser 
based on a Nd:YV 04 gain medium operating at a repetition 
rate stabilized at 100 MHz by an atomic clock and at the 
wavelength of 1064 nm. Each pulse is upconverted with a 
PPLN crystal to a wavelength of 532 nm and energy 1 nJ. 
The pulses, after the MZI, are sent to the Coude path of the 
MLRO telescope, that directs the state |flyout) toward the satel¬ 
lite while actively tracking its orbit. We selected three satel¬ 
lites in low-Earth-orbit (LEO) - Beacon-C, Stella and Ajisai - 
which are equipped with efficient cube-corner retroreflectors 
(CCR). Thanks to the CCR properties, the state is automati¬ 
cally redirected toward the ground station, where it is injected 
into the same MZI used in the uplink. 
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The value of ip{t) originating from the satellite motion can 
be precisely predicted on the base of the sequence of measure¬ 
ments of the instantaneous distance of the satellite, or range 
r, which is realized in parallel. The range is measured by a 
strong Satellite Laser Ranging (SLR) signal at 10 Hz and en¬ 
ergy per pulse of 100 mJ. Thanks to an atomic clock, the SLR 
pulses are separated precisely by AT = 100 ms and synchro¬ 
nized with the 100 MHz pulses. By measuring the tempo¬ 
ral separation AT' of the SLR pulses at the receiver after the 
satellite retroreflection, it is possible to determine the instan¬ 
taneous satellite velocity relative to the ground station Vr{t). 
Indeed, since by the Doppler effect AT' = ^-j^AT, the ve¬ 
locity Vr{t) can be estimated as Vr{t) = C at'+at - 
aration AT' is related to the range r by AT' = AT + Arjc, 
where Ar is the variation of the satellite distance between two 
subsequent SLR pulses. Then, by measuring the range ev¬ 
ery 100 ms, the instantaneous satellite velocity relative to the 
ground station Vr{t) can be estimated, from which (p{t) can 
be derived by Eq. (|^. In the top panel of Fig. for a given 
passage of the Ajisai satellite, we show the measured satel¬ 
lite distance and the estimated ip{t) as a function of time from 
the beginning to the end of the satellite tracking. Since Vr{t) 
is continuously changing along the orbit, the value of (p{t) is 
varying accordingly. In the bottom panel of Fig. 2 we show 
the variation of the theoretical output probability Pc{t) along 
the Ajisai orbit as predicted by eq. 0. 

By the synchronization technique described in Appendix [B] 
we determined of the expected (fret) and the measured (fmeas) 
instant of arrival of each photon. In this way, the histogram 
of the detections in the temporal window of 10 ns between 
two consecutive pulses as a function of the temporal differ¬ 
ence A = tineas — tref Can be obtained. In Fig. [^we show 
such histograms corresponding to constructive and destructive 
interference in the case of satellite Beacon-C. 

In particular, for the constructive interference. Fig. [3^, 
we selected the detections corresponding to (p (mod 27r) € 
[47r/5, 67r/5]. For the destructive interference. Fig. 1^, we 
selected a kinematic phase ip (mod 27r) G [—tt/S, tt/^. The 
detections in the central peak are respectively higher or lower 
than the sum of the two lateral peaks in the two cases. We 
note that the peak width is determined by the detector timing 
jitter which has standard deviation a — 0.5 ns. These two 
histograms clearly show the interference effect in the central 
peak. On the contrary. Fig. 3C is obtained by taking all the 
data without any selection on (p. In this case, the interference 
is completely washed out. These results show that, in order to 
prove the interference effect, it is crucial to correctly predict 
the kinematic phase ip imposed by the satellite motion. 

By using the data of Fig. 3, we experimentally evaluate the 
probability as the ratio of the detections associated the 

central peak N^. to twice the sum Ni of the detections associ¬ 
ated to the side peaks, namely 


p(exp) _ 

" 2Ni ■ 


(3) 


The values = 0.87±0.I0 and = 0.20±0.03 are 

obtained for constructive and destructive interference respec¬ 
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FIG. 3. Constructive and destructive single photon interference 

(Beacon-C satellite, 11.07.2015 h 1.33 CEST). (A) Histogram of sin¬ 
gle photon detections as a function of time A = tmeas — t^ef real¬ 
ized by selecting only the returns characterized by ip{ mod 27r) € 
[Itt/S, Ott/S] that lead to constructive interference. Solid line shows 
the tri-Gaussian fit. By evaluating the Gaussian integrals we ob¬ 
tained the counts Ni = 112 ± 11 for the sum of lateral peaks and 
Nc = 196 ± 14 for the central one. (B) Histogram of single photon 
detections realized by selecting only the returns characterized by 95( 
mod 27r) £ [—tt/S, tt/S]. Here Ne = 112 ±11 and W = 46 ± 7. 
(C) Histogram of single photon detections without any selection on 
the phase. As expected, interference is completed washed out and we 
measured Nc = 1245 ± 35 and Ne = 1306 ± 36, fully compatible 
with Pc = 1/2. In all panels, dotted red lines represent the expected 
counts in case of no interference. 


tively. The values deviates with clear statistical evidence from 
0.5, which is the expected value in the case of no interference. 

A more clear evidence of the role of ip{t) can be demon¬ 
strated by evaluating the experimental probabilities as 

a function of ip. Fig. 4 shows for ten different values of 

the kinematic phase ip and for the three different satellites. By 
fitting the data by = I(I — Vexp cos(/i), we estimated 

the experimental visibilities Vg^p = 67 ± 11% for Beacon-C, 
Vexp = 53 ± 13% for Stella and Vexp = 38 ± 4% for Aji¬ 
sai. The data were collected at the following satellite distance 
ranges: from 1600 to 2500 km (Ajisai, 12.07.2015, h 3.42 
CEST), from 1100 to 1500 km (Stella, 12.07.2015, h 3.08 
CEST) and from 1200 to 1500 km (Beacon-C, 11.07.2015, 
h 1.33 CEST), giving two-way channel lengths ranging from 
2200 up to 5000 km. The interference patterns in Fig. 4 
clearly demonstrate that the coherence between the two tem¬ 
poral modes is preserved along these thousand kilometer scale 
channels with rapidly moving retroreflectors. We attribute the 
different visibilities to residual vibrations of the unbalanced 
MZI between the upgoing and downgoing pulses, since the 
intrinsic visibility of the interferometer was measured to be 
above 95%. To improve the visibility it would be necessary to 
redesign the full interferometric setup to further mitigate this 
effect. We note that, in the double-pass configuration, the in- 
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FIG. 4. Experimental interference pattern. Experimental prob¬ 
abilities as a function of the kinematic phase measured for 

three different satellites. By fitting the data we estimate the visi¬ 
bilities Vexp = 67 ± 11% for Beacon-C, Vexp = 53 ± 13% for 
Stella and Vexp = 38 ± 4% for Ajisai. Dashed lines correspond 
to the theoretical value of Pc predicted by eq. 0. The points 
are obtained by considering ten intervals of the phase defined by 
Xj = [i^Tv, For each interval we selected the data corre¬ 

sponding to (/p( mod 27r) £ Xj: from such data we determined the 
experimental probability of detection in the central peak and 

we averaged the corresponding phase (fi. We note that at point ip = 0 
and (fi = 2 ti the same subset of data were selected. 


terferometer is sensitive to vibrations with frequencies higher 
than 1/rtt (itt is the round trip time). Indeed, the lower visi¬ 
bility is obtained with Ajisai, the satellite with larger distance 
from the ground (from 1600 to 2000 km), with a rtt between 
10.7 ms and 16.7 ms. For the other two satellites the rtt is 
typically lower than 10 ms. 

The mean number of photon /r in the received pulses may be 
derived by measuring the detection rate and using the optical 
losses p ~ 0.27 in the receiving setup. At the primary mirror 
of the receiving telescope, the average ^ during the data acqui¬ 
sition are given by p. « 7 • 10“^ for Beacon-C, p « 2 • 10“^ 
for Ajisai and p « 9 • 10“"^ for Stella. From these values 
we may conclude that interference was probed at the single 
photon level. Indeed, the probability of having more that one 
photon per pulse in the receiver MZI is pp^/2. We estimated 
that that the mean number of photons psat leaving from the 
satellites are Psat*^*'^^ < 20 and < 60 EOll (the in¬ 

stantaneous values fluctuate due to pointing error and turbu¬ 
lence) with a total downlink attenuation between 60 dB and 
70 dB (including the total detection setup losses 20 dB 

due to the optical losses 77 ^ 0.27 and the fact that we used a 
single PMT with 10% efficiency at the output of the MZI). 

Conclusions - Interference at the single photon level be¬ 
tween two temporal modes was observed along a path that in¬ 
cludes a rapidly moving retroreflector on a satellite and with 
length up to 5000 km. We have experimentally demonstrated 


that the relative motion of the satellite with respect to the 
ground induces a varying phase that modulates the interfer¬ 
ence pattern. This varying phase is not present in the case of 
fixed terminals. The effect resulted from the measured inter¬ 
ference pattern during passages of three satellites, Beacon-C, 
Stella and Ajisai, having different relative velocities and dis¬ 
tances from MLRO ground station. 

Up to now, photon polarization was the only degree of free¬ 
dom exploited in long distance free-space Quantum Com¬ 
munications, along ground links of 143 km Il23lj25]l or in 
efforts towards quantum key distribution in space channels, 
of length as large as 2000 km, as recently demonstrated by 
our group El- Indeed, time-bin was never implemented for 
QC over a long-distance free-space channels, fearing that tur¬ 
bulence effects on the wavefront may spoil the interference. 
Here we have demonstrated that atmospheric turbulence is not 
detrimental for time-bin encoding in long distance free-space 
propagation. Indeed, the two temporal modes separated by 
a few nanoseconds are identically distorted by the propaga¬ 
tion in turbulent air, whose dynamics is in the millisecond 
scale I25I : the key point here is the careful matching of the 
interfering wavefronts in the two arms, as shown in Fig. 1 and 
5(A). The results here presented attest the feasibility of time- 
bin/phase encoding technique in the context of Space Quan¬ 
tum Communications. 

Furthermore, the measurement of interference in Space is 
a milestone to investigate one of the big unresolved puzzle 
in Physics, namely the interplay of Quantum Theory with 
Gravitation. As recently proposed theoretically by M. Zych 
et al. ||26l|27l (the optical version of the original Colella- 
Overhauser-Werner (COW) experiment realized with neu¬ 
trons ll28ll l. interference with single photons in Space is a wit¬ 
ness of general relativistic effects: gravitational phase shift 
between a superposition of two photon wavepackets could 
be highlighted in the context of large distance quantum op¬ 
tics experiment. In the case of our setup, the gravitational 
shift for the Ajisai satellite corresponds to about 2 mrad (see 
eq. (23) of lIZTl l. We point out that, unlike the case of ef¬ 
fects manifested by photon polarization rotation, such small 
gravitational effects may get better highlighted by using true 
single photons and by increasing the temporal separation of 
the two interfering modes (with the increase of the difficulty 
in stabilizing the interferometers with large imbalancement). 
To reveal the effect of Gravity in quantum experiments sev¬ 
eral other proposals have been presented lIZTll : these include 
the exchange of elementary particles from moving and ac¬ 
celerated reference frames, which would allow to test Bell’s 
inequalities ll^ [Mil and wavefunction collapse and possible 
gravity-induced decoherence liSTl in laser interferometry with 
a long baseline. 

The interference patterns measured in the present experi¬ 
ment demonstrate that a coherent superposition between two 
temporal modes holds in the photon propagation and its inter¬ 
ference can be indeed observed over very long channels in¬ 
volving moving terminals at fast relative velocity. We believe 
that the results here described attest the viability of the use of 
temporal modes of light for fundamental tests of Physics and 
Quantum Communications around the planet and beyond. 
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Appendix A: Experimental design 

We give a detailed description of the unbalanced Mach-Zehnder interferometer used for the generation and measurement of 
the interference. The interferometer, realized by two beam splitters (BSs), is schematically shown in Fig. 1^. The difference 
between the short (S) and long (L) arm is approximately 1 m. 

In order to observe interference it is necessary to perfectly match the wavefronts of the wave packets traveling through the 
short or the long arm. This matching is particularly necessary for the photon beams reflected by the satellite that are subjected to 
the atmosphere distortion. To this purpose, we exploited two 4/-systems that realize an optical relay of length equal to the path 
difference between the long and short arm of the MZI. Each 4/-system is composed by two lenses with focal length / = 125 
cm, as shown in Fig. |^. It is worth noticing that a single 4/-system, while realizing an optical relay of length 4/, realizes also 
a mirror transformation on the wavefront. A second 4/-system is thus necessary to compensate this mirror transformation. 

In Fig. 1^ we show the performance of the two 4/-systems. At the output of the MZI we imaged the primary mirror while the 
telescope was pointing at a bright star. The two images in Fig. are obtained by blocking the long and short arm respectively. 
The two images are comparable, showing that the wavefronts traveling the short or long arm are well matched at the output of 
the interferometer. 



FIG. 5. (A) Scheme of the unbalanced Mach-Zehnder interferometer used in the experiment. (B) Images of the primary mirror with only the 
short or long arm opened. We exploited the secondary mirror spider and the condensation spots on the primary mirror surface for checking the 
alignment of the interferometer. 


Appendix B: Temporal synchronization 

The quantum measurement at the receiver requires a very precise temporal synchronization and a strong rejection of the 
background, necessary to observe the single photon interference. The measurement of the satellite position along the orbit by a 
laser ranging technique as in Fig. 2 of the main text was also essential for the determination of the expected instant of arrival 
(fret) of the photons at the MZI output. By using a time-to-digital converter with 81 ps resolution (QuTAU), we acquired the start 
and stop signals of the laser ranging pulses, together with PMT detections. We denote PMT timestamps as frneas- In this way, 
we may calculate the histogram of the returns in the temporal window of 10 ns between two consecutive pulses as a function of 
the temporal difference A = fmeas — fret corresponding to a desired value of 


Appendix C: Interference by a mirror moving at constant velocity. 


Here we report the detailed calculations used to obtain equations (1) and (2) of the main text. Indeed, we will derive the effect 
of the reflection by a mirror moving at constant velocity on the superposition between two photon temporal modes. 

Let’s consider a single photon wavepacket whose peak passes at r = 0 for f = 0. Its wavefunction is given by 


where we have defined the quantities 





(Al) 


^ , 

T± = - ± f , 

c 


(A2) 
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and Wo = ^ is the angular frequency. In our convention r is the direction of propagation. The parameter Tc represents the 
coherence time of the pulse. We have used the convention of m , namely 


Tc 


IsMPdr, 


(A3) 


r+co 


9{t) = + t)) = 


+ r)df = e 


(A4) 


If the pulse passes through the unbalanced Mach-Zehnder interferometer shown in Fig. |^, the output is given by 


Mt-) = ^[Mt-) - Mt- - At)] = 


t1 

- e~" 


(A5) 


where At is the MZI unbalancement. 

We now consider a mirror on a satellite moving at constant velocity with respect to the interferometer that at time t = 0 is 
located at r = rgat- We may change reference frame by setting the origin at the location of the satellite. The corresponding 
Lorentz transformation are given by 


r' = 7 (r - Tsat - Pet) 
c 


r = Tsat +^{r' + Pet') 
t = r^it' +P-) 


and where P has been defined after Eq. 3 and 7 = (1 — /3^)“2 is the Lorentz factor. In the upward path, the state ^/>i(r_) is 
directed by the telescope toward the mirror and the wave function in the mirror reference frame can be derived by the transfor¬ 
mations of the T± parameters: 


a\ ! I /1 A /3 / t’sat 

r±= 7 (l±fflr* + —= W—r± + — 


(A 6 ) 


In the mirror reference frame, the mirror reflection can be simply described by the transformation t'_ —t'j^- We now use 

( |A 6 | l, namely ('r+ — to go back to the interferometer reference frame. If we define f^tt = the total 

transformation can be summarized by 


boost to mirror ref. frame /-, r>\ / . ^sat reflection 

-^ 7(1 - P)r_ +- 

c 


- 7(1 - /3)t; + ^ 
c 


boost back to ground ref. frame 


//3('^+ ^rtt) 


where 

It is worth noticing that also normalization should be changed in order to preserve normalization. The beam coming back from 
the satellite is then written as 


V’ 2 ('r+) = V^V'i(-//3('r+ -irtt)) 

7(1 - P) 


V2 


[^o(-//3('r+ - frtt)) - '|/'o(-//3(t+ - Ltt) - Af)]. 


(A7) 


We now comment the two terms. The first term is 


V'o(-//3(t+ - Ltt)) = (/ -^e 


2 -ttLI 






(A 8 ) 


representing a pulse with a coherence time stretched (or compressed) by the Doppler effect: its coherence time is indeed 


, _ Tc _1+ P 
fp ~ I-P^^- 


(A9) 
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FIG. 6. Spacetimes diagrams of light propagation. Short and Long arm pulses are represented by blue and red color respectively. Dashed 
green line represents the satellite trajectory. As explained in the text, pulses separated on ground by a delay At are received with a delay 
At' — Atj fp due to the motion of the satellite. The round trip time (rtt) is Uw Satellite distance r and time t are not to scale. 


We also note that the parameter trtt represents the time that the pulse peak takes to come back to the origin (namely the round 
trip time). The second term in ( |A7| l represents the pulse of eq. ( |A8| ) delayed by At' = The same relation applies also to 
the SLR pulses, that are separated in time by AT = 100 ms. As detailed in the main text, by measuring the temporal separation 
AT' of the SLR pulses at the receiver, it is possible to determine the satellite velocity. Spacetime diagrams for the pulses going 
upward to and downward from the satellite are shown in Fig. |^for different parts of the satellite orbit. 

After passing again in the interferometer we get at the detection port of the MZI (see Fig. 1 of the main text) the following 
state: 


1p3{T+ + trtt) 


»7(1 - /3) 

2 


V^o(-//3T+) + V’o(-//3(t+ + At)) 


- '0o(-Af - fpT+) - - fi3{T+ + At)) 


We used the convention that a BS introduces a i phase shift for the reflected beam and does not change the transmitted beam. 
We now have three pulses at the detector: the probability of getting the photon in the central pulse at r = 0 is given by 


with 


and 


P,{t) = I dt'\M-Mt' + At)) - - fpt')f^ 


111- I dt'^e 


- [1 - V(t) cosifit)] 


+ ^ At) 


2 iuJo{l-f^)At 


ifiit) = a;o[l - f0)]At = , 


V{t) = J- dt'e 


+ _ (t' + At)^ 


= exp{—27r 


At l3{t) 

Tc l + /3(t) 


}■ 


(AlO) 


(All) 


(A12) 
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The above results give equations (1-2) of the main text. 
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